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0\ ■ Abstract 

In the unsplittable flow problem on a path, we are given a capacitated path P and 
n tasks, each task having a demand, a profit, and start and end vertices. The goal is to 
' compute a maximum profit set of tasks, such that for each edge e of P, the total demand of 

Q I selected tasks that use e does not exceed the capacity of e. This is a well-studied problem 

C/3 ' that has been studied under alternative names, such as resource allocation, bandwidth 

O I allocation, resource constrained scheduling, temporal knapsack and interval packing. 

We present a polynomial time constant-factor approximation algorithm for this prob- 
! lem. This improves on the previous best known approximation ratio of O(logn). The 

^ ' approximation ratio of our algorithm is 7 4- e for any e > 0. 

£^ I We introduce several novel algorithmic techniques, which might be of independent 

' interest: a framework which reduces the problem to instances with a bounded range 

I of capacities, and a new geometrically inspired dynamic program which solves a special 



case of the maximum weight independent set of rectangles problem to optimality. In the 
^-^ I setting of resource augmentation, wherein the capacities can be slightly violated, we give 

a (2 -|- e)-approximation algorithm. In addition, we show that the problem is strongly 
NP-hard even if all edge capacities are equal and all demands are either 1, 2, or 3. 



! 1 Introduction 



In the Unsplittable Flow Problem on a Path (UFPP), we are given a path P = {V,E) with 
an integral capacity Ug for each edge e £ E. In addition, we are given a set of n tasks T 
where each task i G T is characterized by a start vertex Si G V, an end vertex ti £ V, a 
demand di G N, and a profit i«j G N. A task i uses an edge e £ E \i e lies on the path from 
Si to tj. The aim is to compute a set of tasks F T with maximum total profit ^i^zpWi 
such that for each edge, the sum of the demands of all tasks in F that use this edge does not 
exceed its capacity. 

The name of this problem is motivated by an interpretation as a multicommodity flow 
problem, where each task corresponds to a commodity. The term "unsplittable" means that 
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the total amount of flow di from each commodity i has to be routed completely along the 
path from the source Si to the sink ti or not at all. There are several settings and applications 
in which this problem occurs, and several other interpretations of the problem. Therefore, 
this problem, and close variants thereof, have been studied under the names bandwidth al- 
location [ini EH [32], admission control [36j, interval packing [22] temporal knapsack [TT] . 
multicommodity demand flow [20], unsplittable flow problem [H [71 [T71 [19], scheduling with 
fixed start and end times [4], and resource allocation [Sj [TU [23[ [36] . In many applications, the 
vertices correspond to time points, and tasks have fixed start and end times. Within this time 
interval they consume a given amount of a common resource, of which the available amount 
varies over time. 

UFPP is easily seen to be (weakly) NP-hard, since it contains the Knapsack problem 
as a special case (in case the path is just a single edge). In addition, Darmann et al. [23] 
show that the special case where all profits and all capacities are uniform is also weakly NP- 
hard. Chrobak et al. [22] strengthen this result by showing strong NP-hardness in this case. 
In addition, they show that the case where the profits equal the demands is strongly NP- 
hard. These results show that the problem admits no polynomial time approximation scheme 
(PTAS) unless P = NP. On the other hand, the special case of a single edge (Knapsack) 
admits an FPTAS. When the number of edges is bounded by a constant, UFPP admits a 
PTAS since it is a special case of Multi-Dimensional Knapsack [27j . 

Most of the research on UFPP has focused on two restricted cases: firstly, the special case 
in which all capacities are equal has been well-studied, which is also known as the Resource 
Allocation Problem (RAP) [SJ [Ml [23l [24l [36] . A more general special case of UFPP is given by 
the No-Bottleneck Assumption (NBA): in that case it is required that maxj di < miug Ug (this 
holds in particular for RAP). We will denote this restriction of the problem by UFPP-NBA. 
For UFPP-NBA, a (2 -|- e)-approximation algorithm is known [20], which matches the earlier 
best approximation ratio for RAP [14| . 

Many previous papers on UFPP partition the tasks into small and large tasks, and use 
different algorithmic techniques for these two groups. For a task i, denote by b{i) the minimum 
capacity among all edges used by task i. For 6 with < 5 < 1, we say that a task i is 6-small if 
di < 6-b{i) holds, and 5-large otherwise. The two main algorithmic techniques that have been 
used in previous results are dynamic programming (for large tasks) and rounding of solutions 
to the linear programming relaxation of the problem (for small tasks) . These techniques work 
well when the NBA holds. 

However, there are several important obstacles that prevent these techniques to be gen- 
eralized to the general case of UFPP. For example, Chakrabarti et al. [T7j show that under 
the NBA the natural LP-relaxation of UFPP has a constant integrality gap. However, with- 
out this assumption the integrality gap can be as large as 0(n) [T7]. Moreover, the NBA 
implies that if all tasks are (5-large, then in any solution there can be at most 2 [l/<5^J tasks 
which use each edge. This property is useful for setting up a dynamic program; see [T7] and 
Section I3.2.2[ Without the NBA this is no longer possible. 

Despite these obstacles, there are a few breakthrough results for (general) UFPP: The 
best known polynomial time algorithm by Bansal et al. [7] achieves an approximation factor 
of O(logn), thus beating the integrality gap of the natural LP-relaxation. This result has 
been generalized to trees by Chekuri et al. [19]. In addition, they gave a linear program- 
ming relaxation for UFPP with integrality gap O(log^re) [19]. Finally, Bansal et al. [6] gave 
a (1 -|- e)-approximation algorithm with quasi-polynomial running time, which additionally 
requires that the capacities and the demands are quasi-polynomial, i.e. bounded by 2^°'^^°^". 
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Nevertheless, it remained an open question whether UFPP admits a constant factor approx- 
imation algorithm (this was asked e.g. in [71 119j). 

1.1 Our Contribution and Outline 

We present the first polynomial time constant-factor approximation algorithm for the general 
case of UFPP. The algorithm has an approximation ratio of 7 -|- e, for arbitrary e > 0. 

To obtain this result we introduce several new algorithmic techniques, which are interesting 
in their own right. We develop a useful viewpoint which allows us to reduce the problem to 
a special case of the maximum weight independent set of rectangles problem. In addition, 
we design a framework which reduces the problem to solving instances where essentially the 
edge capacities are within a constant factor of each other. The techniques can be applied and 
combined in various ways. For instance, for practical purposes, we also show how our results 
can be used to obtain a constant factor approximation algorithm with a reasonable running 
time of only 0(n^ logn). We now go into more detail about these results, the new techniques 
we introduce, and give an outline of the paper. 

Similar to many previous papers, for our main algorithm we partition the tasks into 'small' 
and 'large' tasks. For the small tasks our main result is as follows: For any 7 > and e > 0, 
we present a (3 -|- e)-approximation algorithm for UFPP in Section [3l for the case where each 
task is (1 — 7)-small. We remark that a similar result was given by Chekuri et al. [19], who 
gave an O (log (l/7)/7'^)- approximation algorithm if each task is (1 — 7)-small. Their result 
also applies to trees. To prove our (3 -|- e)-approximation, we introduce a novel framework in 
which the tasks are first grouped into smaller sets, according to their b{i) values, such that the 
techniques for the NBA case can be applied. So the resulting sets can be solved via relatively 
standard dynamic programming, LP-rounding, and network flow techniques. (This is similar 
to e.g. [141 120].) Solutions to these smaller sets leave a small amount of the capacity of each 
edge unused. In our framework we recombine these solutions into a feasible solution for all 
tasks. 

Using the techniques developed for the (3-|-e)-approximation, in Section [3^ we also give a 
result for UFPP in the setting of resource augmentation. We give an algorithm that computes 
a (2 -|- e)-approximative solution which is feasible if we increase the capacity of each edge by 
a factor of 1 -|- for arbitrarily small e > and /3 > 0. Note that this algorithm works with 
arbitrary task sets and does not require the tasks to be small. 

For our main approximation algorithm, it remains to handle the large tasks. For these, 
we present the following main result: for any integer A; > 2, if all tasks are ^-large, we 
give a 2A;-approximation algorithm in Section [H This is based on a geometric viewpoint 
of the problem: we represent UFPP instances by drawing a curve in the plane determined 
by the edge capacities, and representing tasks by axis-parallel rectangles, that are drawn 
as high as possible under this curve. The demand of a task determines the height of its 
rectangle, and the profit of a task determines the weight of the rectangle. Using a novel 
geometrically inspired dynamic program, we show that in polynomial time, a maximum weight 
set of pairwise non-intersecting rectangles can be found. Such a set corresponds to a feasible 
UFPP solution. In addition, we show that when every task is ^-large, this solution yields a 2k- 
approximative solution for UFPP. With this dynamic program we contribute towards the well- 
studied problem of finding a Maximum Weight Independent Set of Rectangles (MWISR) [H 
[30l I34|. Below we discuss this problem in more detail. 

For our main result, we partition the tasks into ^-small tasks and ^-large tasks. For the 
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first group, we apply the aforementioned (3 + e)-approximation algoritfim. For tlie second 
group, our second algoritlim gives a 4-approximation. Returning the best solution of the two 
yields the (7 + e)-approximation algorithm. The main algorithm is summarized in Section [3 
In addition, in Section[5]we show how our results can be combined to obtain a 0(n^ log n) time 
constant factor approximation algorithm for UFPP, and we discuss how our results carry over 
to the generalization from a path to cycle networks, where we obtain a (8 + e)-approximation 
algorithm. 

Finally, we give an alternative proof of the strong NP-hardness of UFPP, which shows that 
a different restriction also remains strongly NP-hard. In the existing NP-hardness proofs \22\ 
[23] , arbitrarily large demands are used in the reductions. In Section [6l we prove that the 
problem is strongly NP-hard even for the restricted case where all demands are chosen from 
{1,2,3} and capacities are uniform (RAP). Note that in contrast to our hardness result, it is 
known that in the slightly more restricted case where the capacities and demands are uniform, 
the problem admits a polynomial time algorithm: In that case, Arkin and Silverberg [3] have 
shown that the problem can be solved in time 0{n^ log n) by minimum-cost flow computations. 

We end in Section [7] with a discussion. In Section [21 notation and terminology are intro- 
duced. First, in the next subsection, we give more background on the many variants of UFPP 
that have been studied in the literature. 

1.2 Related Results 

As mentioned above, the restrictions of UFPP where demands are uniform (RAP) and 
where the No-Bottleneck Assumption maxj dj < miug Ue holds (UFPP-NBA) have been well- 
studied, with the current best approximation algorithm for both problems being the (2 + e)- 
approximation of Chekuri et al. [20]. Both RAP and UFPP-NBA have been generalized in 
various ways: in a scheduling context, one may allow more freedom for choosing the start time 
Si and end time ti of a given task i. Philips et al. [36] obtain a 6-approximation algorithm for 
such a generalization of RAP, by using LP-rounding techniques. For a similar generalization, 
where for each task one out of a set of alternatives needs to be selected, Bar-Noy et al. [8] 
provide a constant factor approximation algorithm using the local ratio technique. RAP and 
UFPP-NBA can be generalized in a different way, more related to network flows, by con- 
sidering graphs other than a path. In graphs other than trees, there may be different paths 
possible between a terminal pair Sj, ti. However, a single path has to be chosen for each se- 
lected terminal pair. This generalization of UFPP to general graphs is called the Unsplittable 
Flow Problem (UFP), or UFP-NBA if the NBA applies. Baveja and Srinivasan [12] provide 
an O ( I i^l )-approximation algorithm for UFP-NBA (on all graphs), improving on various 
earlier results. A simpler combinatorial algorithm with the same guarantee was subsequently 
given by Azar and Regev [5]. Chakrabarti et al. [17] also give an approximation algorithm 
for all graphs. In addition, they observed that an a-approximation algorithm for UFPP-NBA 
gives an (a + 1 + e)-approximation for UFP-NBA on cycles, for any e > 0. This way they 
gave the first constant factor approximation algorithm for both UFPP-NBA and UFP-NBA 
on cycles. Chekuri et al. [20] obtain a 48-approximation for UFP-NBA on trees. 

In addition, many hardness results are known for UFPP (not necessarily under the NBA) 
generalized to various graphs: for general graphs, it is hard to approximate within a fac- 
tor of rid-El^"*^) unless P = NP [5], and for depth-3 trees the problem is APX-hard [29]. 
Hardness-of-approximation results are known even for the special case with unit demands 
and capacities (the Edge Disjoint Path Problem); see [2 [3]. 
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When viewing UFPP as a packing problem, the corresponding covering problem has also 
been studied [151 [16]. that case, tasks have costs instead of profits, and the objective is 
to find a minimum cost set of tasks F, such that for each edge, the sum of the demands of 
all tasks in F that use this edge is at least its capacity. Recently, Chakaravarthy et al. |15j 
designed a primal-dual 4-approximation algorithm for this problem. 

Recall that we reduced UFPP for large tasks to a special case of the Maximum (Weight) 
Independent Set of Rectangles (M(W)ISR) problem. In this problem, a collection of n axis- 
parallel rectangles is given and the task is to find a maximum (weight) subset of disjoint 
rectangles. For the unweighted case of this problem, a randomized 0(log log n)-approximation 
is known [18]. For the weighted case, there are several 0(log n)-approximation algorithms [H 
[30l I34j . The algorithm by Erlebach et al. |26j gives a PTAS for the case that the ratio between 
height and width of rectangles is bounded (note that this does not apply in the special case 
that we need here for approximating UFPP). Our new dynamic programming technique might 
be useful for further research on this problem. 

We remark that this approach for the large tasks is closely related to another well-studied 
variant of RAP: In adjacent resource scheduling problems, one wants to schedule a job on sev- 
eral machines in parallel which must be contiguous, that is, adjacent to each other. In other 
words, this is a variant of MWISR where rectangles are allowed to move vertically, within a 
given range. Duin and van Sluis [21] prove the decision variant of scheduling tasks on contigu- 
ous machines to be strongly NP-complete. RAP on contiguous machines has been considered 
under the name storage allocation problem (SAP), in which tasks are axis-aligned rectangles 
that are only allowed to move vertically. Leonardi et al. [32] provide a 12-approximation 
algorithm for SAP and Bar- Yehuda et al. [9] present a (2 + e + l/(e — l))-approximation 
algorithm. 

In this paper, we study UFPP in the setting of resource augmentation. This means that we 
find a solution which is feasible if we increase the capacity of each edge by a modest factor of 
(l+e). The paradigm of resource augmentation is very popular in real-time scheduling. There, 
the augmented resource is the speed of the machines. For instance, it is known that the natural 
earliest deadline first policy (EDF) is guaranteed to work on m machines with speed 2 — 1 /m 
if the instance can be feasibly scheduled on m machines with unit speed [35]. In addition, a 
matching feasibility test is known [13] . For further examples of resource augmentation results 
in real-time scheduling see [251 EI] • 

2 Preliminaries 

We assume that the vertices of the path P = {V,E) are numbered V = {0, . . . ,m}, and 
E = {{i, i + 1} I < i < m — 1}. We assume that the tasks are numbered T = {1, . . . , n}. 
Recall that tasks are characterized by two vertices Si and tj with Si < ti, and positive integer 
demand di and profit Wi. 

For each task i £ T we denote hy Pi E the edge set of the subpath of P from Sj to tj. 
If e G Pi, then task i is said to use e. For each edge e we denote by Te C T the set of tasks 
which use e. For a set of tasks F we define its profit by w{F) := '^i^pWi. Our objective 
is to find a set of tasks F with maximum profit such that X^jgirnTe — each edge 

e. For each task i we define its bottleneck capacity b(i) by b{i) := miuggp. Ug. An edge e is 
called a bottleneck edge for the task i if e G -Pj and Ue = b{i). In addition, we define for every 
task i that l{i) := b{i) — di. The value l{i) can be interpreted as the remaining capacity of 
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a bottleneck edge of i when i is selected in a solution. Consider a vertex v £ V and an edge 
e £ E with e = {x,x + 1}. We write v < e {or v > e) if v < x (resp. v > x + 1). For 
two edges e = {x, x + 1} and e' = {x' , x' + 1} in E, we write e < e' if x < x' and e < e' if 
X < x' . In other words, we interpret an edge {x,x + 1} simply as a number between x and 
X + 1. Without loss of generality, we will assume throughout this paper that tig > 1 for all 
edges e and di > 1 for all tasks i; zero demands and capacities can easily be handled in a 
preprocessing step. Moreover, observe that one can easily adjust any given instance to an 
equivalent instance in which each vertex is either a start or an end vertex of a task. Such an 
adjustment can be implemented in linear time and it hence does not dominate the running 
times of the algorithms presented in this paper. Therefore, we will henceforth assume that 
m < 2n. Throughout this paper, we will use the notations defined above to refer to the 
UFPP instance currently under consideration. In the few cases where we consider multiple 
instances, the notations will be clear from the context. 

We define an a-approximation algorithm for a maximization problem to be a polynomial 
time algorithm which computes a feasible solution for a given instance such that its objective 
value is at least ^ times the optimal value. Throughout this paper, for a subset of the tasks 
F C T, OPT{F) denotes an optimal solution for the UFPP instance restricted to the task 
set F. The following simple fact shows how we can combine our approximation algorithms 
for different task subsets into one algorithm for all tasks. 

Fact 2.1 Consider a UFPP instance with task set T, and a partition {Ti,T2} of T . If for 
i = 1,2, there exists an ai- approximation algorithm for the instance restricted to the tasks in 
Ti, then there exists an (ai + 02) -approximation algorithm for the entire instance. 

Proof: For i = 1,2, let ALGi denote the solution returned by the approximation algorithm 
for the instance restricted to the tasks in Tj, and OPTi = OPT[Ti) an optimal solution. Let 
OPT = OPT(T) denote an optimal solution for all tasks. So for i = 1,2, w{ALG{Ti)) > 
^w{OPT{Ti)). The algorithm that returns the maximum profit solution of ALGi and ALG2 
has an approximation ratio of qi + a2, since 

w{OPT) = w{OPT n Ti) + w{OPT n T2) < w{OPTi) + w{OPT2) 

< aiw{ALGi) + a2w{ALG2) 

< aimax{w{ALGi),w{ALG2)} + a2ma.x{w{ALGi),w{ALG2)} 

< [ai + Q2) max{w{ALGi),w{ALG2)}. 

□ 

3 Small Tasks 

In this section we present a (3 + e)-approximation algorithm for any set of tasks which are 
(1 — 7)-small (for arbitrarily small e > and 7 > 0). In our main algorithm (for a general 
set of tasks) we will invoke this algorithm as a subroutine for all tasks which are ^-small. 
Moreover, with a slight adjustment of introduced techniques we construct a polynomial time 
algorithm computing a (2 + e)-approximative solution for the entire instance (not only small 
tasks) which is feasible if the capacities of the edges are increased by a factor 1 + e (resource 
augmentation), see Section [331 
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Our strategy is to define groups of tasks such that the bottleneck capacities of all tasks in 
one group are within a certain range. This allows us to compute a feasible solution for each 
group, whose profit is at most a factor 3 + e' smaller than the profit of an optimal solution for 
the group. In addition, each computed solution leaves a certain amount of capacity of every 
edge unused. We devise a framework which combines solutions for a selection of groups into 
a feasible solution for the entire instance, in a way that yields a (3 + e)-approximation (with 
an appropriate choice of e' for the given e). 

3.1 Framework 

We define the framework sketched above. We group the tasks into sets according to their 
bottleneck capacities. Let £ G N be a constant. We define F*^'^ := {i G T\2'' < h{i) < 2'^+^} 
for each integer k. Note that this includes negative values for k, and that at most i ■ n 
sets are non-empty (only those will be relevant later). In the sequel, we will present an 
algorithm which computes feasible solutions ALG (f'^'^) C F^'^ . These solutions will satisfy 
the following properties. 

Definition 3.1 Consider a set F^'^ and let a > and /3 > 0. A set F C F^'^ is called 
{a, /3)-approximative if 

• w{F) > ^■w{OPT{F'''^)), and 

• SieFnTe ^ "Ue — /3 ■ 2^^' for each edge e such that H F^'^ ^0. (In particular, it is a 
feasible solution.) 

An algorithm which computes (a, /3)-approximative sets in polynomial time is called an {a, (3)- 
approximation algorithm. We call the second condition the modified capacity constraint. 

Our framework consists of a procedure that turns an (a, /3)-approximation algorithm for 

each set F^^^ into a (yj-j^ • -approximation algorithm for all given tasks, where q and (3 are 

chosen such that f3 = 2^^'^ > 0. Later, for our resource augmentation result (see Section [3^3]) 
we will work with (a, 0)-approximative sets and therefore, some of the claims below will be 
proven more generally, allowing /3 to be zero. 

Lemma 3.2 (Framework) Let £ £ N and q £ N be constants and let f3 := 2^~'^. Let the 
sets F^'^ be defined as stated above for an instance of UFPP. Assume we are given an {a,j3)- 
approximation algorithm for each set F^'^ with running time 0{p{n)) for a polynomial p. 

Then there is a ^^y^ • -approximation algorithm with running time 0{m • p{n)) for the set 
of all tasks. 

Now we describe the algorithm that yields Lemma 13.21 Assume that we are given an 
(a, /3)-approximation algorithm which computes solutions ALG (F^'^) C F^'^. The key idea 
is that due to the unused edge-capacities of the sets ALG (f'^'^) , the union of several of these 
sets still yields a feasible solution. With an averaging argument we will show further that 
the indices k for the sets ALG (F^'^) that we want to combine can be chosen such that the 

resulting set is an ^^^^ • -approximation. Formally, for each offset c G {0, ...,£ + q — 1} we 
define an index set i]{c) = {c + i ■ (i + q) \ i £ Z} (the values i and q that it depends on 
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will always be clear from the context). For each c G {0, ...,£ + q — 1} we compute the set 
ALG{c) = IJfcer?{c) ALG{F^'^). We output the set ALG{c*) with maximum profit among all 
sets ALG{c). In Lemma we will prove that the resulting set is an (^^^ • a)-approximation. 
First, in Lemma 13.41 we will prove that each set ALG{c) is feasible (using that /3 > 2^"'^). 
This requires the following property. 

Proposition 3.3 Let F be a feasible UFPP solution such that b[i) < u for all i £ F. Then 
for every edge e, YlieFnT^ < 2n. 

Proof: Let e be an edge. If Ue < 2u, then the claim follows immediately. Now suppose that 
Ue > 2u. Any task i G F CiTf, must use an edge whose capacity is less than u. In particular, 
it must use either the closest edge to the left of e or the closest edge to the right of e 
whose capacity is less than u. The total demand of tasks in F using is less than u, and 
the same holds for e/j. It follows that the total capacity used by tasks in F n Tg must be less 
than 2u. □ 



Lemma 3.4 Let ^ G N and q e N be constants. For each set F^'^ let ALG (f'^'^) be a 
{a, j3)- approximative set with (3 > 2^^'^ > 0. Then for each c G {0, + q — 1} the set 
ALG{c) = Ufe677(c) ALG{F^^^) is feasible. 

Proof: Consider a set ALG{F^'^). By definition of the (a, /3)-approximative sets, ALG{F^'^) 
leaves [3 ■ 2^ > 2^~^^~'^ units of the capacity of every used edge free. Observe that this 
is at least twice the maximum bottleneck capacity of tasks in F^~(^+i)'^, Therefore, by 
Proposition [331 the set ALG{F''^^) U ^LG(F''-(^+3)'^) is feasible. In fact, it again leaves a 
fraction of the capacities free, which makes it possible to continue this argument for further 
sets with i > 2, and prove that their union is feasible. 

Formally, let c G {0, ...,£ + q — 1} and let e be an edge. Denote P := I + q. Let k be the 
largest integer in ry(c) such that 2^ < Ue. For every x G f?(c), denote 

i&Ter\ALG{F^-^) 

Since ALG{F^'^) is an (a, /3)-approximation and 2^""^ < /3, we have that 

Ul <Ue- f3-2^ <ue- 2^+1-^. 

For every i > 1, applying Proposition 13.31 for the modified capacities u'^ = Ue — (i • 2^~'^'^ 
shows that 

Summarizing, we have that 

jeTenALG{c) 



i=0 



k-i-P 
e 



i=l 



< Ue + E 2^'+i-(*-i)-^-^ - E 2^+1-^-^-5 
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□ 



Lemma 3.5 Let £ € N and q € N be constants. For each set F^'^ let ALG (^F'^'^) be a 
(a, (3) -approximative set with /3 > 0. Then for the offset c* which maximizes w{ALG{c)) = 
Efc6»?{c) w{ALG{F^^^)) it holds that w{ALG{c*)) > ■ ^-wiOPT), where OPT denotes an 
optimal solution of the given instance of UFPP. 



Proof: Every task is included in £ different sets F^'^. Using this fact, we calculate that 

i+q-l e+q-1 

w{ALG{c)) > 2^ - • ^ {opt {f^'^ 

c=0 c=0 fc6»)(c) " 

^•5^?i;(0Pr(F^'^)) >'^-Yw(0PTr\F^'') =^-w{OPT). 



So there must be a value c such that w{ALG{c)) > ■ ^ ■w{OPT). In particular, this holds 
for c*. □ 

Proof of Lemma \3.S^ Lemma 13.41 shows that ALG{c*) is feasible and Lemma 13.51 shows 

that ALG{c*) is a {j-j^ • -approximation. For computing ALG{c*) we need to compute 

the set ALG (^F^'^^ for each relevant value k. There are at most mi G 0{m) relevant values 
k. Finding the optimal offset c* can be done in 0{m) steps. This yields an overall running 
time of O (m • p{n)) (recall that p{n) is the polynomial bounding the running time needed to 
compute the sets ALG {F^'^)). □ 



3.2 An Approximation Algorithm for Small Tasks 

Now that we have developed the framework to translate (a, /3)-approximation algorithms for 
the sets F^'^ into an approximation algorithm for the entire instance (Lemma I3.2p . it remains 
to present such an (a, /3)-approximation algorithm. In this section, we present a (2 + 
approximation algorithm for sets F^'^ in which all tasks are (1 — 2/?)-small (for arbitrarily 
small /3 > 0). Together with our framework of Lemma 13.21 this yields a (3 + e)-approximation 
algorithm for UFPP for the case that all tasks are (1 — 7)-small, for arbitrary e > and 7 > 0. 
To get some intuition, the reader can think of /3 being equal to min{e,7/2}. 

Suppose we are given a set F^'^ with only (1 — 2/3)-small tasks. In order to derive the 
mentioned (2 + jr^, /5)-approximation algorithm, we choose a value 6 > and split the set 
into 5-small tasks {tiny tasks) and tasks which are 5-large but (1 — 2/3)-small {medium tasks). 
We define S such that for the tiny tasks there is a (y^, /?)- approximation algorithm, presented 
in the following subsection. For the medium tasks, we give a (2, /3)-approximation algorithm 
in Section [3.2.21 



3.2.1 An Approximation Algorithm for Tiny Tasks 

We show that for given e > and /? > 0, there is a 5 > such that if all tasks are (5-small, 
then for each set F^'^ there is a (y^, /5)- approximation algorithm. The key idea is to use 
linear programming techniques and a result by Chekuri et al. [2D] about the integrality gap 
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of the canonical LP-relaxation of UFPP under the no-bottleneck assumption (NBA). UFPP 
with a task set T can be formulated in a straightforward way as an integer linear program: 

IP : max Wi ■ Xi 

s.t. Xi ■ di < Uf> y e & E 

Xie{0,l} ViGT 

The LP relaxation is obtained by replacing the constraint Xj G {0, 1} by < Xj < 1. Chekuri et 
al. [2D] gave an algorithm for UFPP instances which satisfy the NBA (i.e. maxj di < miug Ue), 
in which all tasks are (^-small, which returns a UFPP solution that is at most a factor f{6) 
worse than the optimum of the LP relaxation. Here f{6) is a function for which the limit is 
1 as 6 approaches zero. In other words: The integrality gap of the canonical LP-relaxation is 
1 + e if the NBA holds and all tasks are sufficiently small. 

This result can be used for sets F'''^: By definition, tasks in F'''^ use only edges e with 
Ue ^ 2*^. Call these relevant edges. It is therefore possible to choose the value 6 small enough 
to ensure that the NBA holds, when considering only the relevant edges and (^-small tasks 
in F^'^. Furthermore, modifying the capacities by choosing u'^, = Ue — (5 ■ 2^ decreases the 
capacities of relevant edges at most by a factor 1 — /?. Therefore, the optimal value of the LP 
relaxation also becomes at most a factor 1 — /3 smaller. These are the key ideas to prove the 
following lemma: 

Lemma 3.6 For every combination of constants e > 0, < /5 < 1, and i £N, there exists a 
S > such that if all tasks are 6-small, then for each set F^'^ there is a (j^, (3) -approximation 
algorithm. 

In the remainder of this subsection we prove Lemma 13. 6i Assume we are given constants 
e > 0, £ G N, and /3 with < /3 < 1. For an instance I of UFPP, we denote by LP{I) the 
natural LP-relaxation of the IP-formulation given above, where each constraint Xi £ {0, 1} 
is replaced by < Xj < 1. By LPopt^^) we denote the optimum value of the LP. We 
define f{5') = ^ ^1/^ s' ' '^^^ following result is proved by Chekuri et al. [20], although an 
exact analysis of the running time is not given. We observe that their algorithm admits a 
O(n'^logn) implementation. 

Lemma 3.7 Consider a UFPP instance I for which the NBA holds, in which all tasks are 
5' -small, with 5' < . Then in time 0{n^logn), a feasible UFPP solution ALG fori can 
be computed with w{ALG) > f{5')~^ ■ LPopril). 



Proof: The algorithm of Chekuri et al. works as follows. The tasks are partitioned into 
at most n groups, depending on their demands. The demands and the capacity are scaled 
such that a problem with uniform demands and uniform capacities is obtained. Since the 
demands and capacities are uniform, this can be solved optimally in time O(n^logn), using 
the algorithm by Arkin and Silverberg [4j. Then Chekuri et al. show that combining the 
solutions of each group yields a feasible solution. Furthermore, they have shown in \20\ 
Corollary 3.4] that the obtained solution is at most a factor f{S') worse than the optimal LP 
solution. □ 
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Lemma 3.8 Let 5 > 0, 5' > 0, /3 > and £ G N &e constants such that 6 < and 
5' = < ■ If all tasks in F^'^ are 5-small, then in time 0{n^\ogn), a solution ALG 

for F^'^ can be computed that is {l^^^ fi)- approximative. 

Proof: Consider the UFPP instance X that is obtained by restricting the instance to the 

tasks in F^'^, and only considering edges used by these tasks. So for every task i, it holds 

that 2^ < b{i) < 2^+^, and di < 6 ■ b{i). For every edge e, we have Ue > 1^ ■ 

Now construct the instance X' from X by modifying the capacities as follows: := 
— /3 ■ 2^^, for each edge e. The instance X' contains the same task set as X, without 

modifications. For a task i, by b{i) and we denote its bottleneck capacity in X and X', 

respectively. Hence = b(i) — /3 ■ 2^. 

We will first argue that the algorithm from Lemma 13.71 may be applied to X'. First, 

we note that since in each task is J-small with respect to the original capacities, under the 

modified capacities each task is still J'-small: 

= • (1 - /3) • m < 5' ■ ib(i) -(3.2'^) = 6'- b'd). 

Recall that we required that 6' < ^—j-^, so this condition of Lemma 13.71 is satisfied. Now we 
argue that for the modified capacities, the NBA holds. Since S < {1 — /3)/2^, for all tasks i 
and edges e it holds that di < 6 ■ 2^'+^ < 2^ - ^ ■ 2'' < Ue - /3 ■ 2^ = u'^. This shows that the 
NBA also holds with respect to the modified capacities. 

Hence we may apply the algorithm from Lemma 13.71 to obtain a feasible solution ALG 
for X', with w{ALG) > f{5')-^LPoPT{X'). We argue that LPopt{X') > (1 - /3) • LPopt{X): 
Since for every edge e it holds that > 2'^ , it follows that u'^ = Ue -13-2^ > {1-I3)ue. Thus, 
if we take a feasible solution to LP{X), and scale all the variables Xi by a factor (1 — /3), we 
obtain a feasible solution to LP{X'), in which the objective value has also been scaled by a 
factor (1 — /3). This gives that 

Lemma [3. 71 , 

w{ALG) > fiS'y' ■ LPopt{X') 



> . (1 - /3) . LPopt{X) 



where OPT{X) is an optimal (integer) UFPP solution for the instance X. □ 
The proof of Lemma 13.61 now easily follows: 

Proof of Lemma \3.(k Observe that the limit of f{6') is 1 as 5' approaches zero. Hence for every 
e>0, 0</3<l and £ G N, we can choose a 5' > and 5 := (1 - (3)5' such that f{5') < 1 + e, 
5 < and 6' < Then Lemma 13.81 shows that if all tasks are 5-small, then in time 



O(n'^logn), a solution ALG for F^'^ can be computed that is (y^, /3)-approximative. □ 



3.2.2 An Approximation Algorithm for Medium Tasks 

It remains now to find a (2, /3)-approximation algorithm for tasks in F^'^ that are both (5-large 
and (1 — 2/3)-small (for the 6 we obtained from Lemma [3.6|) . for arbitrarily small /3 > 0. When 



A Constant Factor Approximation for Unsplittable Flow on Paths 



12 



restricting to sets F^'^ with only 5-large tasks, the essential property is that for any edge e 
and any feasible solution F, there are at most ^(^t) ^^^^s in F that use e. This property 
allows for a straightforward dynamic program to be used to compute an optimal solution 
(see e.g. [HI HZ]). This can be turned into a (2, /3)-approximate solution: since tasks are 
(1 — 2/3)-small, it can be shown that in polynomial time, any solution F can be partitioned 
into two sets which are both feasible for the modified capacities. Using these ideas, we will 
prove the following lemma in the remainder of this subsection. 

Lemma 3.9 Let f3 > 0, 6 > and I £ N be constants and assume we are given an instance 
of UFPP in which all tasks are both (1 — 2j3)-small and 5-large. Then, for each set F^'^ there 
is a {2, 13) -approximation algorithm. 

Suppose we are given a set F^'^ whose tasks are all (5-large and (1 — 2/3)-small. 

Proposition 3.10 Let F C F^'^ be a feasible solution in which all tasks are 5-large. Then 
for any edge e, at most 2^~^^ /5 tasks in F use e. 

Proof: Let F be a feasible solution. For each task i G F^'^ it holds that b{i) > 2^^. In 
addition, all tasks are (5-large. Therefore, for all tasks i £ F it holds that di > 5-2^. For 
every task i £ F, b{i) < 2}^^^ . So according to Proposition 13.31 for every edge e it holds that 
X^ieFnTe ^« ^ 2^^+^+^. Therefore, at most = tasks in F use e. □ 

Due to Proposition 13.10] there is a straight forward dynamic programming algorithm for 
computing an optimal solution [141 117j . (We remark that the running time is vP^ ^^).) 

Proposition 3.11 For constant 5 > and £ G N, if all tasks are 5-large, an optimal solution 
for F^'^ can be computed in polynomial time. 

Proof sketch: By Proposition 13.10] each edge can be used by at most 2^/5 tasks. Hence, for 
each edge e there can be at most n'^(^ combinations of tasks which use e in an optimal 
solution. We enumerate all these combinations for each edge e. For each of these, we establish 
a dynamic programming cell which stores the maximum profit one can obtain from tasks i 
with Si < e, given the respective combination of tasks that use e. The correct values for 
these cells can be computed by iterating through the edges of the path from left to right. 
See [a HZ] for details. □ 

Lemma 3.12 Let /3 > and let F be a feasible solution for F^'^ in which all tasks are 
(1 — 2/3) -small. Then in time 0{in?), F can be partitioned into two sets and which are 
both feasible for the modified capacity constraints u'^, := Ue — P ■ 2^ . 

Proof: Note that the claim is trivially true if /3 = by setting := F and := 0, so 
now assume /3 > 0. We initialize two sets := 0, := 0. Assume that the tasks in F are 
ordered such that the start vertices are non-decreasing. We consider the tasks in F in this 
order. In the i-ih. iteration we take the task i. We add i to a set with I G {1, 2} such that 
U {i} obeys the modified capacity constraint, i.e., it leaves a free capacity of /3 • 2^' in each 
edge. 
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It remains to show that indeed either U {i} or U {i} obeys the modified capacity 
constraint. Assume to the contrary that neither U {i} nor U {i} obey the modified 
capacity constraint. Then there are edges e,e' on the path of i such that 

di>>Ue-(3- 2^ (3.1) 

i'e(//iu{i})nTe 

and 

^ di, >Ue'- p-2^. (3.2) 

i'G{H2u{i})nT,, 

Inequality (|3.ip implies that J^i'&H^nTe^^' ^ Ug — (5 ■ 2^ — di. Inequality (|3.2p gives that 
^i'&H^nT , ^i' > Ue' — (3 ■ 2^^ — di- Assume w.l.o.g. that e < e' or e = e'. Recall that we 
considered the tasks by non-decreasing start index. Hence, all tasks in {H^ U {i}) n Tg' use e 
as well. For the next calculation we need that 

di<{l- 2p)h{i) < Ue'(l - 2/3) 

and hence Uf,' — di > 2/3 ■ Ugi. Also, note that Ug/ > 2'^ since i S F C F^'^. We calculate that 

Ue > 
> 

> 
> 

This is a contradiction. Hence, task i can be added to one of the sets such that U {i} 
still obeys the capacity constraint. When computing and we need to check for each 
task in F whether adding it to one of the sets violates the modified capacity constraint in one 
of the edges. Since w.l.o.g. m < 2n, this check can be done in 0{n) time for each task. There 
are n tasks in total, and hence the entire procedure can be implemented in O(n^). □ 

Now we can prove Lemma 13.91 

Proof of Lemma \3.9\: Since the tasks are 5-large, we can compute an optimal solution 
OPT{F^'^) in polynomial time (Proposition 13. IT]) . Since tasks are (1 — 2/3)-small, in time 
O(re^), the solution OPT{F^'^) can be partitioned into two solutions and that obey 
the modified capacity constraint (Lemma I3.12p . Returning the solution of these two with 
maximum profit then yields a (2, /5)-approximation for F^'^. □ 

3.2.3 An Approximation Algorithm for (1 — 7)-Small Tasks 

By combining Lemmas 13.61 and 13.91 we obtain an (a, /3)-approximation algorithm for small 
tasks. 



^ d^/ + di, + di 

i'eH^nTe / \i'£H2r\T^, / 
Ue -p-2^- di^ + (ue' -P-2^-di)+d^ 

Ue + Ue' - -2^-2'' 
Ue + 2l3Ue' -2(3-2^ 

Ue. 
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Lemma 3.13 Let e > 0, /? > 0, and £ be constants with /3 < 1 and assume we are given 
an instance of UFPP in which all tasks are (1 — 2j3)-small. Then, for each set F^'^ there is 
0^ (2 + (3) -approximation algorithm. 

Proof: Given e, /? and i, Lemma [3. 61 shows that there exists a 5 > such that for the 5-small 
tasks we have an (yz^, /3)-approximation algorithm (for each set F^'^). The remaining tasks 
are both 5-large and (1 — 2/3)-smaU. If 5 > 1 — 2/3, we are done. Otherwise, Lemma 13.91 
shows that we have an (2, /3)-approximation algorithm for these (for each set F^'^). Together 
this gives an (2 + /3)-approximation algorithm for each F^'^ (observe that Fact 12.11 also 
applies to (a, /3)-approximation algorithms). □ 

Using the above lemma and our framework (Lemma 13. 2p . we obtain our main result for 
small tasks. 

Theorem 3.14 For any pair of constants e > and 7 > 0, there is a polynomial time 
(3 + e)- approximation algorithm for UFPP instances in which all tasks are (1 — ^) -small. 

Proof: Choose ^ G N, g e N, /3 := 2^"'? and e' > such that 1-2/3 > 1 -7 (hence all tasks are 
(1 — 2/3)-small), and ■ ^2 + yr^^ < 3 + e, which is always possible^. Now, combining the 

framework using the chosen i and q (Lemma l3.2p with (2 + /3)-approximation algorithms 
for every F^'^ (Lemma 13. 13p yields a (3 + e)-approximation algorithm. □ 

3.3 Resource Augmentation 

Using the techniques derived above, we now describe a polynomial time algorithm that com- 
putes a set of tasks T' C T such that w{T') > (2 + e)~^ • w{OPT), and T' is feasible if the 
capacity of every edge is increased by a factor 1 + e. Note that for this result we do not require 
that the tasks T are (1 — 7)-small. 

The main idea is the following: from the results in Sections 13. 2.1 1 and l3. 2. 21 we will conclude 
that there are (1+e, 0) and (1, 0)-approximation algorithms for the tiny tasks and all remaining 
tasks, respectively. Combining these, we obtain a (2 + e, 0)-approximation algorithm for each 
set F^'^ (without any further conditions on its tasks) in Proposition 13.151 We apply our 
framework with the sets computed by this algorithm. In Lemma [3.16l we show that the union 
ALG{c) = Ufce77(c) ALG{F^'^) is feasible when the capacities of the edges are increased by a 
factor 1 + 2^^'' (this lemma takes the role of Lemma 13.41 from the original framework). 

The first step is to establish the approximation algorithms for the sets ALG{F^'^). 

Proposition 3.15 Let e > 0. There is a (2 + e, ^)- approximation algorithm for each set F^'^ 
which runs in polynomial time. 

Proof: Using Lemma 13.61 with /? = yields that for all e > there is a 5 > such that 
there is a (1 + e, 0)-approximation algorithm for each set F^'^ which consists only of (5-small 

^ For instance, choose j3 small enough such that 1/(1 — /?) < 1 + e', 1 — 2/? > 1 — 7, and that there is an 
integer q > 2 with /3 — 2^^''. Now choose £ £N such that -^i^ < 1 + e'. We obtain an approximation factor of 
(2 + ii^) • ^±2 < (2 + (1 + e')^) • (1 + e') which is at most 3 + e if e' is sufficiently small. 
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tasks. Proposition 13.111 implies that for each fixed 6 > there is a polynomial time (1,0)- 
approximation algorithm (i.e., an optimal algorithm in the usual sense) for sets F'''^ which 
consist only of (5-large tasks. Using Fact 12. ll this yields a (2 + e, 0)-approximation algorithm 
for each set F'''^. □ 

The next lemma is an adjusted version of Lemma l3.4i In contrast to the latter, here the 
increased capacity allows us to combine the computed sets ALG{F^'^) to a globally feasible 
solution, without requiring that every solution should leave a fraction of the capacity free. 
Recall that r/(c) = {c + i ■ {I + q)\i & Z}. 

Lemma 3.16 Letq^n,^e N, and c G {0, + 1}. For each k, let ALG^'^ C F^'^ he a 
feasible solution and define ALG{c) = Ufeer;(c) ALG{F^'^). Then for each edge e it holds that 

d,<ne- (1 + 22-9). 

ieALG{c)nTe 

Proof: Let k be the largest integer in r]{c) such that 2'^ < Ug. Let P := i + q. From 
Proposition 13.31 we conclude that "^i^ALC^'^nTe <2 - 2'^+^ for each k < k. This implies that 



oo 



d^ = Y "^i+Yl Yl 



< Ue + 2 • 2 



p+i 
1=1 

oo 

= Ue + 2^+^+^X^(2^)- 
i=l 



e 2P_1 

= Ue { I -\ t: 

\^ 2^+9 - 1 , 

< Ue- (1 + 22-9). 

□ 

The following lemma establishes our framework (see Lemma [3.2p in the setting of resource 
augmentation. 

Lemma 3.17 Let I be a UFPP instance. Let ^ G N and q £ N be constants. If polynomial 
time [a, 0)- approximation algorithms exist for each set F^'^ , then there is a polynomial time 
algorithm that computes a set of tasks whose profit is at least ■ ■ w{OPT{X)), and 
which is feasible if the capacity of every edge e is modified to ite(l + 2^-9) . 

Proof: The algorithm works as follows. For each c G {0, + q — 1} we define ALG{c) := 
UfcGr;(c) ALG{F^'^). We output the set ALG{c*) with maximum profit among all sets ALG{c). 
Lemma[33]impliesthat i(;(^LG(c*)) > j^^-^-w{OPT{l)). LemmaEHDyields that ALG(c*) 
is feasible if the capacity of each edge e is changed to U(.{1 + 22-9). □ 

The previous lemmas can be combined to obtain the main result of this section: 
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Theorem 3.18 Let e > 0, /3 > 0. Given an instance I of the UFPP problem. There is 
a polynomial time algorithm that computes a set of tasks whose profit is at least (2 + e)"-"^ • 
w{OPT[X)), and which is feasible if the capacity of every edge e is modified to We(l + /?)■ 

Proof: We choose ^ G N, g G N, and e' > such that 2^^^ < /3 and ^ • (2 + e') < 2 + e, 
which is always possible (similar as in the proof of Theorem I3.14p . The claim follows by 
combining the framework of Lemma 13.171 with the (2 + e', 0)-approximation algorithms from 
Proposition I3.15[ □ 



4 Large Tasks 

In this section we provide a polynomial time 2fc-approximation algorithm for instances con- 
sisting of only l/fc-large tasks. In our main algorithm (in Section E]), this is used as a 4- 
approximation algorithm for the set of 1/2-large tasks. The main idea is to restrict to UFPP 
solutions of a certain form: we will represent tasks by rectangles drawn in the plane, and 
compute an independent set of rectangles of maximum weight (profit). 

By (xi, yi, X2, 2/2) we will denote the axis-parallel rectangle in the plane with upper left 
point (xi,?/i) and lower right point {x2,y2)- We will call two rectangles (xi, yi, 2:2, 2/2) and 
(x'^, y'^, Xg, 2/2) compatible if they do not overlap (i.e. do not share an internal point). More 
precisely, they are compatible if at least one of the following holds: X2 < x[, x'2 < xi, yi < 
or y[ < y2. 

Definition 4.1 (associated rectangle) With a taski we associate the rectangle {si,b{i),ti,i{i)). 

Note that Si and ti are integers since the path vertices are labeled 0, . . . , m, and that b{i) 
and i{i) := b{i) — di are non-negative integers as well. Tasks are called compatible if their 
associated rectangles are compatible, and a task set F is called an independent task set (ITS) 
if all tasks are pairwise compatible. 




Figure 4.1: A UFPP instance consisting of a path of length eight (with capacities 
2,4,4,6,6,5,3,3) and four tasks, represented using rectangles drawn as high as possible un- 
der the capacity profile. 

The following geometric interpretation motivates this choice of rectangles, see Figure WA\ 
For every edge e = {x, x -|- 1} € E{P), we draw a horizontal line segment between (x, u^) and 
(x -|- l,ne). We add vertical line segments to complete this into a closed curve from (0,0) to 
(m,0), called the capacity profile (represented as a bold grey curve in Figure \4A\ and later 
figures). For every task i, the associated rectangle is now a rectangle of height di, drawn as 
high as possible under the capacity profile. The tasks in Figure UTT] have profits 3, 2, 2, and 2, 
respectively. Therefore, an optimal UFPP solution consists of tasks 1, 3, and 4, with a total 
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profit of 7. However, since the rectangles for the tasks 1 and 4 are incompatible, this is not 
an ITS. The optimal ITS consists of tasks 2, 3 and 4, with total profit 6. This ITS is a UFPP 
solution as well. This is in fact always true, as shown in the next proposition. 



Proposition 4.2 Let F QT be an ITS. Then F is a feasible UFPP solution. 



Proof: Consider an edge e G E{P), and all tasks Tg that use e. Consider two tasks i, j € FnT^. 
They are compatible, but Si < e < tj and sj < e < ti, so either l{i) > or > h{i) must 
hold. It follows that the tasks in F n Tg can be numbered ii,. . . ,ip such that 6(ii) < £{12), 
^(^2) < ^(^3)) etc. Hence, 



p 



Y,d^, = Y.^h{ik) - ^{ik)) 
k=l k=l 



p-1 

b{ip) - i{ii) +yz{b{ik) - i{ik+i)) 



k-1 



< Up 



□ 

In Section [4. 11 we will first show that if all tasks are ^-large, then the profit of an optimal 
ITS is at most a factor 2k worse than the profit of an optimal UFPP solution. We will prove 
this by showing that any UFPP solution can be partitioned into 2k ITSs. 

After that we will give a dynamic programming algorithm for finding a maximum profit 
ITS. The core concept is that of a corner, which corresponds to a region of the plane that 
contains a subset of the rectangles. In S ection [4 . 2 1 corners are defined, and a recursion formula 
is given for computing the profit of an optimal ITS that fits in such a corner. The dynamic 
programming algorithm that is based on this recursion is subsequently given in Section [4.31 



4.1 Partitioning UFPP Solutions into ITSs 

In this section we prove that if all tasks are 1 //c-large, the value of an optimal ITS is at most 
a factor 2k worse than the value of an optimal UFPP solution, by showing that any UFPP 
solution can be partitioned into 2k ITSs. This result is based on the following property. 



Proposition 4.3 Let F be a feasible UFPP solution, and let k > 2 be an integer, such that 
every task in F is 1/k-large. Let e be a bottleneck edge for i £ F. Then there are at most 
k — 1 tasks j £ F \ {i} that are incompatible with i and use e. 



Proof: Let the set F' <^ F \ {i} consist of all tasks that are incompatible with i and that 
use e. Suppose > k. Consider j G F' . Since j and i are incompatible, but both 
use e and e is a bottleneck edge for i, it follows that b{j) > i{i) = b[i) — di. Therefore, 
di + J2j&F' dj >di + l YjjdF' Kj) > di + ^\F'\(b{i) - di) > b{i) = Ue, contradicting that F is 
a feasible solution. □ 

A partition {Fi, . . . , Fi} of a task set F into I ITSs will be encoded by an £-coloring a of 
F. This is a function a : F ^ {1, . . . ,£} such that if a(i) = a(j), then i and j are compatible. 
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(So a{i) = j means that i G Fj.) An edge e is called a separator for a task set F if it is a 
bottleneck edge for some task i £ F, such that all tasks in F that use e are incompatible 
with i. So by Proposition 14.31 if e is a separator edge, then there are at most k — 1 tasks in 
F that use e and are incompatible with i, when F is a UFPP solution consisting of l/fc-large 
tasks. A coloring a of F is called nice if for every edge e and task i £ F that has e as 
its bottleneck edge, all tasks that use e and are incompatible with i are colored differently. 
The main idea behind these notions and our construction of a coloring is as follows: We will 
identify a separator edge e, and consider the set Fq of tasks i with Sj < e, and Fi of tasks i 
with ti > e. (Note that FqUFi = F and Fq n Fi / 0.) Unless Fq = F or Fi = F, we may 
use induction to conclude that both admit a nice 2A;-coloring. Then, since e is a separator 
edge and the colorings are nice, all tasks in Fq D Fi are colored differently in both colorings. 
Therefore these can be combined into a single nice 2fc-coloring for F. However, since it may 
occur that Fq = F or Fi = F, we need a slightly more sophisticated argument, which we will 
now present in detail. 

Lemma 4.4 Let F be a feasible UFPP solution, and let k >2 be an integer, such that every 
task in F is 1/k-large. Then there exists a nice 2k-coloring for F. 

Proof: We prove the lemma by induction over \F\. The statement is trivially true when 
\F\ < 2k. Now suppose that \F\ > 2k, and assume that the above statement holds for every 
such set F' with < \F\. The proof is illustrated in Figure HT2l Let be a bottleneck 
edge for some task in F, with minimum capacity among all such edges. Let L C F be the set 
of tasks that use cb- By Proposition 14.31 and the choice of es, \L\ < k. 

Let Es = {ei, . . . , Cp} be the set of edges e for which there is a task i with e as bottleneck 
edge, such that i S L or i is incompatible with some j G L. (So cb G Es-) The edges in Es 
are numbered such that if x < y, then Cx < Cy. 

We first argue that all edges in Es are separators for F. For cb £ Es, the statement is clear 
since cb is a bottleneck edge with minimum capacity. Now consider an edge e G Es \ {cb}, 
and let z G F be a task with e as bottleneck, that is incompatible with some task i' G L. 
Then < b{i') holds. Now suppose there is a task i" using e that is compatible with i. 
Since both i and i" use e and e is a bottleneck edge of i, from their compatibility it follows 
that b{i") < £{i) < b{i') = cb- But this contradicts that is a bottleneck edge for F with 
minimum capacity. 

Now, for 1 < j < p — 1, define Fj C F to be all tasks i with Si < Cj+i and ti > Cj. 
Similarly, define Fq <^ F to be all tasks i with Sj < ei, and define Fp <Z F to be all tasks i 
with ti > Cp. Observe that FqU . . .U Fp = F. 
Case 1: For every j, \Fj\ < \F\. 

In this case, we may use induction to conclude that for every Fj there exists a nice 2k- 
coloring aj. We can combine these into a nice 2A;-coloring of F: all tasks in C = FqCi Fi use 
the edge ei. Since ei is a separator, there exists a task i £ C with ei as bottleneck such that 
all tasks using ei are incompatible with i. So, in both the nice 2A;-coloring ao of Fq and the 
nice 2A;-coloring ai of -Fi, all tasks in C are colored differently. Therefore, we can permute 
the colors of ai such that the tasks in C receive the same colors in ao and ai. At this point, 
the colorings can be combined into a 2A;-coloring a' of FqU Fi, which is again nice. 

Next, we can combine the nice 2fc-coloring 02 of F2 with the nice 2A;-coloring a' of FqUFi in 
a similar way, and continue like this with F3, . . . , Fp, until a nice 2A;-coloring of F = FqU. . .UFp 
is obtained. This proves the desired statement in this case. 
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Case 2: There exists a j such that Fj = F . 

For such a choice of j, define F'- = Fj \L. Since L is non-empty, \Fj\ < \F\ holds, so we 
may again use induction to conclude that admits a nice 2A;-coloring a. We will extend 
this to a nice 2A;-coloring of Fj = F. For ease of exposition we will assume that 1 < j < p 
(the cases j = and j = p are similar but easier). So both ej and e^+i are edges in Es- We 
assume also w.l.o.g. that e_B > ^j+i- 

We observe that every task i G that is incompatible with some task in L uses either ej 
or ej+i: by definition of Es, such a task i must use some edge G Es- If x > j + 1 then by 
definition of Fj, i also uses Cj+i. If x < j then similarly i also uses ej. 

Since cb > Cj+i and Fj = F, all tasks in L use Cj+i. Then, since e,+i is a separator, 



there are at most k — \L\ tasks in Fj that use ej+i (Proposition 14. 3| ): let Fr denote these 
tasks. In addition, there are at most k tasks in Fj that use ej; denote these by F^. Since 
\Fl\ + |Fr| < 2k — \L\ and there are 2k available colors, we can choose \L\ different colors 
that are not used for tasks in F^ U Fr, in the coloring a. We extend the nice 2A:-coloring a 
for Fj to Fj = F by using these \L\ colors for the tasks in L. Since we observed that all tasks 
in Fj that are incompatible with some task in L are included in Fl U Fr, it follows that this 
yields again a nice 2A;-coloring for F. Since we now constructed a nice 2A;-coloring for F in 
all cases, this concludes the proof of Lemma l4.4i □ 




Figure 4.2: An illustration of the proof of Lemma 14.41 for the case k = 2 (using logarithmic 
scale on the y-axis, all tasks are ^-large). The gray task is the single task in L, and the bold 
tasks are associated with the indicated separator edges ei, . . . , 65. 

We will now show that this bound is tight for every k > 2: consider a path P of length 
5, with capacities 2k'^, 2k'^ + 2k, 2{2k'^ + 2k), 2k'^ + 2k and 2/c^, in order along the path. (So 
V{P) = {0, . . . , 5}.) Introduce k — 1 tasks with Sj = 0, = 3 and di = 2k + 1, k — 1 tasks 
with Si = 2, ti = 5 and di = 2k + 1, one task with Sj = 1, = 3 and di = 2k + 3, and one 
task with Si = 2, ti = A and di = 2k + 3. All tasks can be verified to be ^-large. In addition, 
they all satisfy b{i) > 2k'^ > £{i) and Sj < 2 < 3 < tj, so all are pairwise incompatible, and 
therefore at least 2k colors are needed. Finally, they constitute a feasible UFPP solution: 
for the first and fifth edge this is easy to see. For the second and fourth edge, the sum of 
demands using the edge is {k-l){2k + l) + {2k + 3) = 2k'^ + k + 2 < 2k'^ + 2k = 2} = -"{3,4}, 
since k > 2. All tasks use the third edge, so the demand sum is 2- {{k — l)(2A;-|-l)-|-(2A;-|-3)) = 
2 • (2A;2 -\- k + 2) = Ak"^ + 2A; + 4 < Ak'^ + Ak = ^{2,3}, since k>2. The above is summarized 
in the following proposition. 
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Proposition 4.5 For any k £N with k >2 there is an instance of UFPP with the properties 
that 

• the instance consists of 2k tasks which are all ^-large, 

• the set of all tasks of the instance yields a feasible UFPP solution, and 

• the optimal ITS contains only one task. 

4.2 Computing Optimal ITSs in Corners 

To bound the complexity of our algorithm, it is useful to assume that all edge capacities are 
different. In the next lemma we first show that by making small perturbations to the capac- 
ities, this property can easily be guaranteed, without changing the set of feasible solutions. 

Lemma 4.6 Let I be a UFPP instance with task set T = In linear time, the 

capacities and demands can be modified such that for the resulting instance I' : 

• All edge capacities are distinct, and 

• a set of tasks F QT is feasible for I if and only if F is feasible for I' . 

Proof: Recall that the m path edges are labeled {i, i + 1}, for i G {0, . . . ,m — 1}. For every 
edge e = {i,i + 1}, we change the capacity to u'^ = mue + i- For every task j G T, we change 
the demand to d'j = mdj. This gives the instance I', in which all capacities are distinct. 

Suppose F C T is feasible for I'. Then for every edge e = + 1}, it holds that 
T^iennrdi = ikYlieT^npd'i < [^u'^\ = [ue + ±\ = Ue- Hence F is also feasible for the 
original instance I. Clearly, every task set F that is feasible for I remains feasible for I'. 
□ 

The central concept of our dynamic program is that of a corner {x,y,z). First we give 
an informal, geometric explanation of this notion, using the representation of the problem 
explained before. Subsequently, we will give a formal definition. 

A corner {x, y, z) corresponds to a certain region under the capacity profile, see Figure 
In our dynamic program, we will compute for every such corner the profit of an optimal 
top-drawn solution that fits entirely within this region. This will be done using previously 
computed values for 'smaller' corners. A corner {x',y\z') is sma/Zer than the corner {x,y,z) 
if its region is a strict subset of the region associated with (x, y, z) (we will make this more 
precise in the proof of Lemma I4.16P . 

A corner is determined by an integer x-coordinate < x < m (a path vertex), and two 
integer y-coordinates y > and z > 0. If the capacity of the edge {x — 1, x} G E{P) is more 
than y, then draw a horizontal line segment from the point (x, y) to the left, to the first point 
that lies on the capacity profile curve. This point will be called {wL,y). Otherwise, if the 
capacity is at most y, then let wl = x. Similarly, if the capacity of {x,x -|- 1} G E{P) is 
more than then draw a horizontal line segment from {x,z) to the right, to the first point 
that lies on the capacity profile. This point will be called {wr,z). Otherwise, let wn = x. 
Connect these line segments into a single curve from {wL,y) to {wii,z) by adding a vertical 
line segment from {x,y) to {x,z). This curve and the capacity profile together now enclose a 
bounded region, which is shown as a shaded region in Figure 14.31 (In the special case that 
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Figure 4.3: The corner C{x,y,z) contains all tasks for which the associated rectangles are 
fully contained in the shaded region. 

> z > u^x~i.x} > y or u^x-~i,x} > y ^ ^{z,z+i} > Zj the corner actually corresponds 
to two disjoint regions.) This is the region shown that we associate with the corner {x,y,z); 
we say that a task fits into the corner (x, y, z) when its rectangle, drawn as explained earlier, 
lies fully in (the closure of) this region. Now we give formal definitions. 

Definition 4.7 A triple {x, y, z) of integers is called a corner if < x < m, y > and z > 0. 
For a corner {x,y,z), we denote by WL{x,y) or simply wl the lowest numbered path vertex 
such that for all i with wl < i < x, it holds that Ujj > y. Similarly, wr{x,z) or simply 
wr is defined to be the largest numbered path vertex such that for all i with wr > i > x, it 
holds that j} > z. Hence, wl < x and wr > x. 

For each corner (x, y, z) we define a set C{x, y, z) which — intuitively speaking — consists of 
all tasks for which the associated rectangles are contained in the region associated with this 
corner. 

Definition 4.8 For a corner {x,y,z), we denote by C{x,y,z) the set of all tasks i for 
which at least one of the following holds: 

• WL{x,y) < Si, ti < wr{x,z) and i{i) > ma,x{y,z}, 

• WL{x,y) < Si, ti < X and £{i) > y, or 

• X < Si, ti < wr{x, z) and i{i) > z. 

Due to this definition, we say that task i fits into the corner {x, y, z) or corner (x, y, z) 
contains i \i i ^ C{x,y,z). Hence, C{x,y,z) is the set of all tasks that fit into the corner. 
Note that they are possibly incompatible. 

For a given UFPP instance and corner {x,y,z), we denote by P{x,y,z) the maximum 
value of w{F) over all ITSs F with F C C{x,y,z). An ITS F with F C C{x,y,z) and 
w{F) = P{x,y,z) is said to determine P{x,y,z). Let Umax := ^^^eeE{P) ^e- Observe that 
all tasks fit into the corner {m,0,Umax), so computing P(m,0, Umax) yields the desired value: 

Proposition 4.9 P{m, 0, Umax) equals the profit of an optimal ITS, where m is the path 
length. 



A Constant Factor Approximation for Unsplittable Flow on Paths 



22 



We will now show how P{x, y, z) can be computed in various cases. Proposition 14.101 
collects the easy cases (which are easily understood using the above geometric interpretation) , 
and Lemma 14.121 covers the complex case which gives the main recursion formula to compute 
an optimal ITS. In Proposition 14.10] statement (i) and (ii) show how to rewrite a corner in 
'standard form' without changing the corresponding region. Statement (iii) deals with the 
case where the corner corresponds to two disconnected regions. 

Proposition 4.10 

(i) Ify = z, thenC{x,y,z) = C{wr{x, z),y,Umi,x), and hence P{x,y, z) = P{wr{x, z),y,Uma.x). 

(ii) If X = or y > u^x_i x}, then C{x,y,z) = C{x,Umax, z), and hence P{x,y,z) = 

P{Xj liinax; -2) . 

(Hi) If y < z, X > I, y < and n|^_i^^} < z < then C{x,y,z) = 

C(x, y, Umax) U C{x, limax, z) and hence P{x, y, z) = P{x, y, Umax) + P{x, 

Proof: The statements (i) and (ii) follow immediately from the definitions. We will now prove 
statement (iii). Consider an ITS F that determines P{x,y,z). Since u^x-i,x} ^ z, every task 
i £ F with Si < X — 1 and ti> x has h{i) < z. Hence F contains no tasks i with Si < x — 1 
and > X + 1, and thus can be partitioned into two ITSs that fit into the corners (x, y, limax) 
and (x, Umax, -2) respectively. It follows that P{x,y,z) < P{x,y,Uma.^) + P{x,Uma.x, z). 

Now let Fi and F2 be ITSs that determine P{x,y,Umi^x) and P(x, Umax, -z)- These are 
disjoint and both fit in the corner {x,y,z), so P{x,y,z) > w{Fi) + w{F2) = P{x,y,Uraax) + 

P{x,Uma.x,z). □ 

Because of Proposition 14.101 (i). we only need to consider corners {x,y,z) where y < z or 
y > z holds. Since these cases are symmetric, we will only consider corners (x, y, z) for which 
y < z holds, in the following definitions and lemmas. One can easily deduce the corresponding 
statements for the case y > z. For the sake of brevity and readability, we will however not 
explicitly treat this case. By Proposition 14. 10) (ii). we may furthermore assume that x > 1 and 
y < u^x-i,x}- Informally, Proposition 14.101 (iii) shows that we only need to consider corners 
for which the associated region is connected. These observations show that we may restrict 
our attention to special types of corners, which we call proper. 

Definition 4.11 A corner {x,y,z) is called proper if y < z, x > 1, y < us^^_i ,^^ and either 
z < U{^_i^^} or z> U{^^^+i} holds. 

Now we will consider the more complex case where the given corner {x,y,z) is proper. The 
main idea is as follows. Consider an ITS F that determines P{x, y, z). Either F also fits into 
the smaller corner (x — 1, y, z), or there exists a task j with < z and tj = x. In the latter 
case, we show that F can be partitioned into two task sets that fit into smaller corners, and one 
single task i. Given a task i with ti < x, we will consider the two smaller corners {si,y,b{i)) 
and {x,b{i),z). These are illustrated in Figure l474l For the indicated task i and the depicted 
top-drawn set F, we have the desired property that C C{si,y,b{i)) U C{x,b{i), z). We 

will show that there always exists a task i for which this holds; such a task is called special. 
This explains the following recursive formula for proper corners. 
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WL{x,b(i))f 

Wfi{si,b(i)) 




{x,z) wn{x,z) 
{x,b{i)) 



WLix,y) 



WLix,y) 



Figure 4.4: On the left, a proper corner {x,y,z) and an ITS F C C{x,y,z) is shown. When 
choosing a special task i £ F, it holds that F \ {i} C C{si,y, b{i)) U C{x, b{i),z). 



Lemma 4.12 Consider a UFPP instance in which all edge capacities are distinct. Let 
{x, y, z) be a proper corner. Then 



P{x,y , z) = max <^ P{x — l,y, z), max 

i£C{x,y,z),ti<x 



^Wi + P{si,y, b{i)) + P{x, b{i), z) } I . 



Before we can prove Lemma [4.121 we will define which properties are required for a special 
task, and prove that such a task always exists. The essential property of a special task i is that 
the rectangle {si,b{i), x,y) is compatible with any associated rectangle of a task in F \ {i}, 
and does not overlap with the capacity profile. In this case we will say that there are no tasks 
in F that lie to the right of i or lie below i, which is defined as follows. 

Definition 4.13 Consider a proper corner {x,y,z), and a task i G C{x,y,z) with ti < x. 

• A task j G C{x,y, z) lies to the right of i if ti < Sj < x and £{j) < b{i). 

• A task j G C{x,y,z) lies below i if b{j) < i{i) and tj > Sj. 



Definition 4.14 Let {x,y,z) be a proper comer and F C C{x,y,z) be an ITS. A task i £ F 
with ti < X is called special with respect to F and (x, y, z) if 

• for all edges e with Si < e < x, it holds that Ue > b{i), and 

• there is no task j £ F \ {i} that lies to the right of i or that lies below i. 

We now sketch how to find a special task i in case that F ^ C(x — l,y,z). We start 
with a task i such that ti = x and £{i) < z holds. Now, assume there is a task j whose 
associated rectangle is incompatible with the rectangle {si,b{i),x,y). Then, j lies below i 
(that is, b{j) < i{i)), and no associated rectangle of a task crosses the line segment from 
{sj,b{j)) to (x,6(j)). This follows essentially from the fact that all tasks are compatible in 
F, and all associated rectangles touch the capacity profile. We continue the same procedure 
with j in the role of i until we find a special task i. We now prove this formally. 

Lemma 4.15 Let {x,y,z) be a proper corner and F C C{x,y,z) be an ITS. Then F C 
C(a; — 1, y, z) or there exists a special task i £ F with respect to F and (x, y, z). 
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Proof: A task i E C{x,y,z) is called a candidate for {x,y,z) if 

• ti < X, 

• i{i) < z and 

• for all edges e G E{P) with Si < e < x, Uf. > b{i) holds. 

From u^^_i,i.j > y it follows that WL{x,y) = wl{x — l,y). Obviously, wr{x — 1,z) < 
'Wr{x,z). Therefore, since y < 2; it holds that C{x — l,y,z) C C{x,y,z). Furthermore, 
wr{.x — 1,2) < wr{x,z) can only occur if < z. But since the corner is proper, this 

implies that u^^ ,j,^iy < z, so wr{x, z) = x. These observations show that the only case when 
F Q C{x — l,y,z) does not hold is when there is a task j £ F with tj = x and < z. 
Hence it is easily seen that j is a candidate task (for {x,y,z)). In addition, since tj = x, no 
task in F lies to the right of j. So to prove the lemma statement, we may assume that there 
exists at least one candidate task in F with no tasks to the right of it. Choose i G -F to be 
such a task with minimum value for b{i). We prove that no task in F lies below i, which will 
prove the lemma. Suppose to the contrary that a task j & F lies below i, so b{j) < < b{i) 
and tj > Si. It is easily checked that j is a candidate as well. By choice of i, it must then 
hold that there exists a task k £ F that lies to the right of j (otherwise j should have been 
chosen in the role of i). So tj < Sk < x, and £{k) < b{j) < i{i) < b{i). Therefore, Sk < ti] 
otherwise k would also lie to the right of i. But now we can use the fact that i and k are 
compatible: Recall that we have that Sk < ti, Sj < tj < Sk, and i{k) < b{j) < b{i). So 
the only case in which i and k can be compatible is when b{k) < < b{i). This however 
contradicts the fact that all edges between Si and x have capacity at least 6(i); the edge that 
determines b{k) lies between Si and x. We conclude that i is a candidate task that satisfies 
the additional properties that there are no tasks in F below it, or to the right of it, which 
therefore is special. □ 

Now we can prove Lemma 14.121 



Proof of Lemma \4-12 : Let {x,y,z) be a proper corner in a UFPP instance in which all 



capacities are distinct. We will first show that 

P{x,y,z) < max \P{x-l,y,z), max \wi + P(si,y,b{i))+P(x,b{i),z)\\. (4.1) 

1^ i&C{x,y,z),ti<x (. J J 

Let F be an ITS that determines P{x,y,z). If F C C(x — l,y,z) then P{x,y,z) < P{x — 
l,y,z). If not, then by Lemma [4.15l there exists a special task i. We argue that -F\{2} can be 
partitioned into two ITSs Fi and F2 with Fi C C (sj, y, b{i)) and F2 C C (x, b{i), z) , which will 

prove that P{x, y, z) < + P(sj, y, b{i)) +P(x, b{i), z) |, and therefore the Inequality (j4.ip . 

First we consider the wl and wr vertices for both corners. Let cb = {v,v-\-l} be a bottleneck 
edge for i, which is unique since all capacities are distinct. So Ueg = b{i) and < cb < ti. 
For the corner {si,y,b{i)), it is easy to see that WL{si,y) = WL{x,y) and Wfi{si,b{i)) = v. 
Now consider the corner {x,b{i), z). Obviously the wr value is the same as for C{x,y,z). In 
addition, since i is special, all edges e with Si < cb < e < x have Ue > b{i). Therefore, since 
we assumed all edge capacities are distinct, Ue > b{i) holds for all such edges e. It follows 
that wl{x, b{i)) = v + 1. 

Now consider a task j £ F \ We distinguish five cases. 
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1. tj < Si: WLix,y) = WL{si,y) and £{j) > y imply that j e C{si,y,b{i)). 

2. Si < tj < cb'- Since j is compatible with i but does not lie below i, > b{i) holds. 
Therefore j G C(sj, y, b{i)). 

3. Sj > x: From j G C{x,y,z) it easily follows that j G C(x,6(i),z). 

4. cb < Sj < x: Since j is compatible with i but does not lie to the right of i, > b{i) 
must hold. In addition, if tj > x then > z. Therefore j G C(x, z). 

5. < cb and < tj-. We argue that this is not possible. We have that < Ueg = b{i). 
Because i and j are compatible, it then must be that b{j) < But this contradicts 
that there is no task below i. 

Since we considered all possibilities, this concludes the proof of Inequality (|4.ip . 
Next, we argue that 



Recall that since {x,y,z) is a proper corner, C{x — l,y,z) C C{x,y,z), so P{x,y,z) > 
P{x — l,y,z). Now consider a task i G C{x,y,z) with ti < x. Let Fi and F2 be the ITSs 
that determine P{si,y,b{i)) and P{x,b{i), z), respectively. We will argue that Fi n F2 = 
and FiL) F2 C C{x,y,z). Since clearly i ^ Fi and i i<2, it will follow that P{x,y,z) > 
Wi + P{si,y,b{i)) + P{x,b{i), z), proving Inequality ()4.2p . 

Now let cb a bottleneck edge for i. Then wji{si,b{i)) < cb- On the other hand, since 
ti < X, it holds that WL{x,b{i)) > cb- It follows that Fi F2 = 0. Finally, we prove that 
both sets fit in the corner {x,y,z). Recall that WR{si,b{i)) < cb < ti < x. Since in addition 
6(i) > y, it follows that C{si,y,b{i)) C C{x,y,z). (Because Wfi{si,b{i)) < x, it is irrelevant 
whether b{i) > z or not.) It is obvious that C{x,b{i), z) C C(x,y,z), since b{i) > y. This 
concludes the proof of Inequality (14. 2p . □ 

4.3 A Dynamic Programming Algorithm for Finding Optimal ITSs 

In this section we give an O(n^) dynamic programming algorithm for computing an opti- 
mal ITS (Theorem I4.18p . The main recursive routine V{x,y,z) for computing P(x,y,z) is 
shown in Algorithm [H which is based on the recursive formulas from the previous section 
(Proposition 14.10] and Lemma l4.12p . At the end of this section, we summarize how this al- 
gorithm, together with Lemma l4.4| yields the 2A;-approximation algorithm for UFPP in the 
case of l/fc-large tasks (Theorem 14. 20p . 

First, we show that the recursion terminates, by showing that all recursive calls are made 
with arguments {x',y',z') that yield a 'smaller' corner than the original argument {x,y,z). 

Lemma 4.16 Algorithm[l\ terminates. 

Proof: We show that when given an input (x, y, z), all recursive calls will be with arguments 
{x',y',z') such that the corner {x',y',z') is strictly smaller than {x,y,z), with respect to the 
following size measure: For a corner (x, y, z) with y < Umax and z < Umax, we define 




max 

i&C{x,y,z),ti<x 




area(x, y, z) = {x - WL){uma.x - v) + {wr - x)(umax - z). 
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Algorithm 1: V{x,y,z): a recursive algoritlim for computing P{x,y,z) 
Input: A UFPP instance where itmax is the maximum edge capacity, all edge capacities 

are distinct, and integers 0<x<m,0<y< Umax) < z < Umax- 
Output: P{x,y,z). 

1 if X = or y > U{^_i^,j.y then y := UmiL^; 

2 ii X = m or z > Uj^^^+i} then z := itmaxj 

3 Compute wl ■= WL{x,y) and wr := wrIxjz); 

4 if Wl = wr then return 0] 

5 li y = z then z := n^ax; x := wr] 

6 if y < z then 

if < z < then return {V{x, yi "Umax ) +P(x,nmax,^;)); 

return max — 1, y, z), 

maxigc(x,y,^),i,<x (^^i +7^(si,2/,&(«)) + 7^(3:, 2:)) }; 
end 

(In the remaining case, y > z holds:); 
9 if < y < U{2.„i .j,} then return (7^(x, y, Umax) + 7^(x, n^ax, ^;)); 

10 return max |'P(x + 1, y, z), 

C(x,?;,z),s,>x +'P{ti, b{i),z) +V{x, y, 6(i)) 



max, 



Note that for all corners (x, y, z) that may be considered, area(x, y, z) is a non-negative integer. 
Hence the lemma statement then follows. 

A recursive call is made in Line [7] whenever a corner (x,y,z) is considered that satisfies 
the following bounds. Firstly y < z < Umax, and therefore x > 1 and y < it{a;_i^a;}. Hence 
Wl < X. If the if-condition is satisfied then wr > x holds. This shows that for both calls 
y, Umax) and V{x,Umax, z), the area strictly decreases. Now consider Line [51 which is 
reached whenever the corner is a proper corner. In this case, when considering the corner 
{x — 1, y, z), the wl value does not change, and the wr value can only decrease. Since y < z 
it then follows that area(x — l,y,z) < aicea{x,y, z). Now consider a task i £ C(x,y, z) with 
ti < X. We have that WR{si,b{i)) < x and b{i) > y, so area(sj, y, 6(i)) < area(a;, y, z). Clearly, 
aTea{x,b{i), z) < area{x,y, z). So in all recursive calls in Line [HI the area decreases. By 
symmetry, the same holds for Lines [9| and [TOl □ 



Lemma 4.17 When given a corner {x,y,z), Algorithmic returns P{x,y,z). 



Proof: We prove the statement by induction over the total number of recursive calls (which 
is possible since Lemma 14.161 shows that the algorithm terminates). That is, we will prove 
that a call V{x,y,z) returns the value P{x,y,z), assuming that this holds for recursive calls 
V{x' ,y' , z') that may be made. 

The modification in Line [1] yields an equivalent corner by Proposition I4.10( ii) . By sym- 
metry, the same holds for Line[2j If wl = wr then C{x,y,z) = 0, so the correct answer is 
returned in Line[ll The modification in Line [5] is correct as well (Proposition 14. lOl fi)). So we 
may assume now that wl < wr, and y ^ z. 
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Consider the case that y < z. This imphes that y < Umax, so x > 1 and y < u^x-i,x}- 
If u^x-i,x} ^ ^ z < u^x,x+i}y then it is a non-proper corner, so Proposition 14.101 (iii) 
shows that Line [7] returns the correct answer. So if Line [8] is reached, the considered corner 
is proper. By Lemma |4.12|, Line [8] returns the correct answer. 

The case that y > z is analog; by symmetric arguments, the answers given in Lines [9] 
and [To] are correct. This shows that Algorithm [1] returns the correct answer in every case. 
□ 

Lemma 14.171 and Proposition 14.91 show that the profit of an optimal ITS can be computed 
by calling Algorithm [T] with the argument (m, 0, Umax)- However, to attain the desired com- 
plexity of O(n^), it is necessary to use a dynamic programming table that stores previously 
computed values, to avoid repeating the same recursive calls. The complexity then follows 
from the fact that this table needs at most entries. 

Theorem 4.18 There is an O(n^) algorithm for computing a maximum profit ITS for a 
UFPP instance. 

Proof: The algorithm is as follows. First, in time 0{n + m), we transform the given instance 
into an equivalent instance, in which all edge capacities are distinct (Lemma 14. 6p . As men- 
tioned in Section [21 in time 0{n + m) we can transform the instance to an equivalent instance 
in which all vertices occur as the start or end vertex of some task, so we may furthermore 
assume that m < 2n. This preprocessing does not change n and cannot increase m. On 
the resulting instance we call Algorithm IT] with, arguments (x, z) — O5 ^max 

), to obtain 

the total profit of an optimal ITS (Lemma I4.17( Proposition 14. 9p . With a straightforward 
extension, we can compute an optimal ITS itself. 

A small modification is necessary to obtain a time complexity of 0{nm'^): whenever any 
value P{x, y, z) is computed during the course of computation, this value is stored in a table. 
Whenever a recursive call to V{x,y,z) would be made, the algorithm instead returns the 
stored value P{x, y, z), if it has been computed before (in a different recursion branch). This, 
together with Lemma 14.161 ensures that for every corner (x, y, z) the routine V{x, y, z) is 
called at most once during the course of computation. (Lemma 14.161 ensures that a given 
corner is not considered more than once in the same recursion branch.) 

We call a corner (x, y, z) relevant if y = or ?/ is equal to the capacity of some edge, and 
if the same holds for z. Observe that when calling Algorithm [1] with a relevant corner as 
argument, then every possible recursive call is with a relevant corner as well. So the total 
number of recursive calls is bounded by the total number of relevant corners. There are at 
most m + 1 possibilities for x, and at most m + \ for y and z (since these need to correspond to 
actual capacities), so the total number of recursive calls is bounded by 0{im?). Now consider 
the complexity of a single call. The computations in Line [3] can be done in time 0{m). 
Ignoring recursive calls. Lines [7| and [9| take constant time. Lines [8l and 1101 take time 0(n); 
for every task i, in constant time one can decide whether i € C{x,y,z) and ti < x (resp. 
Si > x), and evaluate the given expressions. The remaining lines clearly take constant time, 
so it follows that a single call of Algorithm [1] takes time 0(n) + 0{m) C 0{n). Hence, the 
total complexity becomes 0{nm^) C O(n^). □ 

Our results can be formulated as follows in terms of the Maximum Weight Independent 
Set of Rectangles (MWISR) problem: 
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Theorem 4.19 MWISR can be solved in time O(n^) for instances with n axis-parallel rect- 
angles, in which each rectangle contains a point {x,y) such that no point {x,y') with y' > y 
is part of a rectangle. 

We summarize the results of Section S] in the following theorem. 

Theorem 4.20 For every integer k >2, there is an O(n^) time 2k -approximation algorithm 
for UFPP restricted to instances where every task is ^-large. 

Proof: In time O(n^) we compute an optimal ITS F for the instance (Theorem I4.18p . This 
is a feasible UFPP solution (Proposition I4.2[ ). Let F* be an optimal UFPP solution. Then 
F* can be partitioned into at most 2k ITSs (Lemma 14. 4p . so w{F) > j^w{F*). □ 

5 Main Approximation Algorithms 

We now apply the main results from the previous two sections to obtain our main approxi- 
mation algorithm. 

Theorem 5.1 For every e > there is a (7 + e)- approximation algorithm for the UFPP 
problem. 

Proof: We partition the tasks into ^-small and ^-large tasks. For the |-small tasks, we have a 
(3 + e)-approximation algorithm (Theorem 13. 14p . and for the ^-large tasks, a 4-approximation 
algorithm (Theorem I4.20p . Returning the best solution of the two then yields a (7 + e)- 
approximation algorithm (Fact 12. ip . □ 

The running time of the above algorithm is dominated by the dynamic program which 
is needed for handling the 'medium size' tasks, see Proposition 13.111 Unfortunately, the 
exponent of the running time is quite large. However, in the following theorem we show 
that already a moderate running time of 0(n^ log n) suffices to obtain a constant factor 
approximation, by partitioning the tasks differently in order to avoid this expensive dynamic 
program. Hence, at the cost of a higher approximation factor, this algorithm has a much 
better running time. 

Theorem 5.2 There exists a 25. 12 -approximation algorithm for UFPP with running time 
0{n^ logn). 

Proof: We partition the tasks into |-small and ^-large tasks. Using Theorem l4.201 we obtain 
an 18-approximation algorithm for the ^-large tasks, that runs in time 0{n'^). For the ^-small 
tasks, we apply our framework with constants ^ = 3, g = 5. According to Lemma 13. 2^ we 
then choose /3 := 2^^'^ = 1/16. We first verify that these choices satisfy the conditions of 
Lemma 13. 8t 

• 5 < (1 - P)/2^ holds, since (1 - /3)/2^ = 15/128 > 15/135 = 1/9 = 6. 

• 6/{l - /3) < (3 - V5)/2 holds, since 6/{l -13) = 16/135 w 0.12 < 0.38 w (3 - y/b)/2. 
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Independent set UFPP instance: 
instance G: 



Left edges 




corresponds to: Ci 62 63 64 65 
capacity: 7676767676 



Right edges 

I'l V2 V3 Vi 

7 7 6 6 4 4 1 1 



vertex: 
Path P: 
Tasks T: 

demand 



12 



2m + 2n 



V2 

demand — 2 



Jl{ 



demand — 3 \ 



demand 



: odd edge 
: even edge 

long high-profit 
short high-profit 
low profit 



Figure 6.1: An example of a graph G with coloring a and the resulting instance ufpp(G). 
(The demand of a task is indicated by the height of its bar.) 



Hence we may apply Lemma 13.81 to conclude that for every integer A:, in time O(n^logn), 
an {/{§')/{! — /S), /?)-approximative solution can be computed for the ^-small tasks in F^'^^ 
where 6' = (5/(1 - /3) = 16/135, and 

nS') = . < 2.503. 



0.5372 



Therefore this is an (a, /3)-approximative solution, with a = 2.503/(1 — f3) < 2.67. Our 
framework (Lemma 13. 2p then gives an approximation algorithm for the instance consisting of 
all g-small tasks, with complexity 0(n^ log n) and approximation ratio • a < 7.12. 

We apply the 18-approximation algorithm for the |-large tasks and the 7.12-approximation 
algorithm for the ^-small tasks. Returning the best solution of these two gives a 25.12- 
approximation algorithm for the entire instance (Fact 12. 1|) . □ 



Chakrabarti et al. [17] showed that any a-approximation algorithm for UFPP implies a 
(1 + a + e)-approximation algorithm for UFP on cycles (also called ring networks). Note that 
on a cycle, for each task i two possible paths can be chosen from Si to tj. (Even though in |17j 
only UFP-NBA was considered, we observe that their argument applies to UFP as well.) This 
yields the following corollary. 



Corollary 5.3 For every e > there is a (8+ e)- approximation algorithm for the Unsplittable 
Flow Problem on Cycles. 



6 Strong NP-Hardness 

In this section we prove that UFPP is strongly NP-hard for instances with demands in {1, 2, 3}, 
using a reduction from Maximum Independent Set in Cubic Graphs. Let G, k be an indepen- 
dent set instance, where G is a connected graph of maximum degree 3 with G ^ K4. The 
question is whether G contains an independent set of size at least k. This problem is NP- 
hard |28j (even for cubic graphs, in which all degrees are exactly 3). Let V{G) = {vi, . . . , Vn}, 
and E{G) = {ei, . . . ,6^}. By Brooks' Theorem (see e.g. [33]), G is 3-colorable since G con- 
tains no K4^. Such a coloring can be found in polynomial time [33], so for our polynomial 
transformation we may assume that a proper 3-coloring a : V{G) —t- {1,2,3} is given. 
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We now construct an equivalent instance of UFPP as follows, see Figure 16.11 The path 
P has 2n + 2m + 1 vertices, labeled 0, . . . , 2n + 2m. For the proofs below it will be useful to 
distinguish different types of edges: an edge {i — € E{P) (with 1 < i < 2n + 2m) is a 
left edge if i < 2m and a right edge otherwise. It is an odd edge if i is odd, and an even edge 
otherwise. Even left edges {2k — 1,2k} (with 1 < A; < m) will be associated with edges 
of G, and even right edges {2k — 1, 2k} (with m + l<k<n + m) will be associated with 
vertices Vk-m of G. 

For every vertex Vi G V{G), introduce the following tasks: First, introduce one long task 
with start vertex and end vertex 2m + 2i — 1. Next, let di, . . . ,ad be the indices of the 
edges incident with Vi, in increasing order. Introduce d+1 short tasks with the following start 
and end vertex pairs: (0, 2ai — 1), (2(Ti, 2a2 — 1), ... , {2ad-i,2ad — 1), (2<Trf, 2m + 2i). For all 
of these tasks for vertex Vi, the demand is equal to a{vi), and the profit is equal to a{vi)n 
times the number of odd edges used by the task. The aforementioned tasks will be called the 
high-profit tasks (for Vi). Finally, for Vi we introduce one low-profit task from 2m + 2i — 1 to 
2m + 2i with profit 1 and demand a{vi). 

Doing this for all vertices gives the tasks of the UFPP instance. It remains to set the edge 
capacities of P: 

• For odd left edges e = {2k — 2,2k — 1} G E{P) for integer 1 < k < m, set Ue = 
Et,ey(G)«(^)- 

• For even left edges e = {2k — 1,2A;} G E{P) for integer 1 < k < m, set Ue = 

(Et,ey{G)«W) - 1- 

• For right edges e = {2k - 2,2k - 1} G E{P) and e = {2k - l,2k} G E{P) for integer 
m + 1 < /c < n + m, set Ue = Yl'i=k-m. 

This gives the instance ufpp(G) of UFPP, to which we will refer in the remainder of this 
section. With T we will denote the set of all tasks of ufpp(G). With an independent set / of 
G we associate a solution of ufpp(G) of the following form. 

Definition 6.1 A set F <^T of tasks o/ufpp(G) is in standard form if: 

1. For every i £ {1, ... ,n} : 

• Either F contains the long high-profit task for Vi and the low-profit task for Vi , but 
no short high-profit tasks for Vi, 

• or F contains all short high-profit tasks for Vi, hut neither the long high-profit task 
for Vi nor the low-profit task for Vi . 

2. For every {vi,Vj} G E[G), F does not contain both the long high-profit task for Vi and 
the long high-profit task for Vj . 

Given an independent set / of G, we can construct a solution F to ufpp(G) by selecting 
the long high-profit task and the low-profit task for Vi if Vi G /, and all short high-profit 
tasks for Vi otherwise. Since / is an independent set, this set F is in standard form; it 
satisfies Property [2] in the above definition. We will first show that the resulting task set F 
is a feasible solution to ufpp(G), by showing that in fact any task set in standard form is a 
feasible solution. Secondly, we show that the task set F constructed this way is an optimal 
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solution for UFPP(G). To this end, we first argue that the total profit of a sohition F in 
standard form depends only on the number of long high-profit tasks in F. Subsequently, 
we will show that any optimal solution for ufpp(G) is in standard form. Together this will 
show that the instances G and ufpp(G) are equivalent (when considering them as decision 
problems with appropriately chosen target objective values). 

Proposition 6.2 Let F T be in standard form. Then F is a feasible solution for UFPP(G). 

Proof: Since F is in standard form, for every Vi £ V{G), every edge of P is used by at 
most one task for Vi. Therefore, for odd left edges e = {2k — 2,2k — 1} G E{P) (with 
1 < k < m), the capacity Ug = J2veV{G) ^(''^) exceeded. Similarly, for right edges 

e = {2k — 2,2k — 1} E E{P) (with m + 1 < k < n + m), the capacity J27=k-m Q^(^i) 
exceeded. Now consider even left edges e = {2k — 1,2k} E E{P) with 1 < k < m. These 
correspond to edges E E{G). Let = {vi,Vj}. Since F is in standard form, it cannot 
contain high-profit tasks for both Vi and Vj. Note that no short high-profit tasks for Vi and 
Vj use the edge {2k — 1, 2k}. Hence the total capacity of e used by F is at most 

max I ^ a{w)^ — a{vi), ^ a{w)^ — a(t>j)| < ^ a{w)^ — 1 = Ue. 

w£V{G) weViG) wevic) 

□ 



Proposition 6.3 Let F Q T be a solution in standard form for UFPP(G), which contains x 
long high-profit tasks. Then w{F) = x -\- Y17=i Oi{vi)n{m -\- i). 

Proof: Recall that for high-profit tasks for a vertex Vi, the profit equals a{vi)n times the 
number of odd edges used by the task. The long-high profit task for Vi (from to 2m -\-2i—l) 
uses all odd edges {2k — 2,2k — 1} of P with 1 < k < m + i. Hence if selected, it contributes 
a{vi)n{m + i) to the total profit. Together, the low-profit tasks for Vi use exactly the same 
set of odd edges of P. So either way, the high-profit tasks for Vi contribute a{vi)n{m + i) to 
the total profit. The single low-profit task for Vi, with a profit of 1, is selected if and only if 
the long high-profit task is selected. This yields the stated total profit. □ 

Lemma 6.4 If G has an independent set I with \I\ = x, then ufpp(G) admits a solution of 
profit X + Yll=i (y-{vi)n{m + i). 

Proof: Construct F by selecting the long high-profit task and the low-profit task for Vi 
if Vi E /, and all short high-profit tasks for Vi otherwise. Since / is an independent set, 
this set F is in standard form. So by Proposition 16. 2| it is a feasible solution to ufpp(G). 
Proposition 16.31 now gives the total profit. □ 

Now we will show that any optimal solution to ufpp(G) is in standard form, which will 
yield the converse of Lemma |6.4[ We say that a solution F fully uses the capacity of an edge 
e E E{P) if J2ieFnTe di = u^. The majority of the total profit of a solution is proportional 
to the capacity use of all odd edges. Since there is at least one solution that fully uses the 
capacity of all of these edges (a solution in standard form), the next proposition follows. 
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Proposition 6.5 Let F be an optimal solution for UFPP(G). For every odd edge e G E{P), 
F fully uses the capacity of e. 



Proof: Recall that for every high-profit task with a demand of d, the profit equals dnx, where 
X is the number of odd edges used by the task. Hence the total profit of high-profit tasks 
in a feasible solution F is bounded by n times the sum of capacities of all such edges. This 
capacity sum is 



4=1 k=l i=k i=l 

Suppose that for at least one odd edge the full capacity is not used by F. The low-profit 
tasks in F can in total only yield a profit of n, so then the total profit is at most 

, n X n 

n ( a{vi){m + i) — 1 I + n = n > a{vi){m + i). 



It >. ft 

a{vi){m + i) — 1 I + n = n^"^ 

1=1 ^ i=\ 



Since there exists a feasible solution for UFPP(G) with strictly higher profit (Lemma 16. 4|) . this 
shows that F is not optimal. □ 

Now we will apply Proposition 16.51 to prove that every optimal solution for ufpp(G) is in 
standard form. 



Proposition 6.6 Let F he an optimal solution for UFPP(G). Then for every i G {1, . . . ,n}, 
F either contains no short high profit tasks for Vi, or contains all short high profit tasks for 

Vi- 



Proof: Suppose that the statement is not true. Then we may assume that there exists an 
even left edge e = {x, x + 1} in P such that F contains the short high profit task for Vi that 
ends at x, but not the short high profit task for Vi that starts at x + 1 (the argument in the 
case where this is reversed is the same). Let ViVj € E{G) be the edge that e corresponds to, 
so Vi and vj are the only vertices of G for which e is not included in the set of corresponding 
short high profit tasks. Since F fully uses the capacity of both of the odd edges {x — 1, x} and 
{x + 1, X + 2} next to e (Proposition 16. 5( ). F must contain tasks starting on x or x + 1, with 
a total capacity of exactly a{vi), to compensate for the fact that the short high profit task 
for Vi starting at x + 1 is not included. But the only such task is a short high profit task for 
Vj, which has a capacity of a{vj) ^ a{vi) (since a is a proper coloring of G), a contradiction. 
□ 



Lemma 6.7 Every optimal solution F for UFPP(G) is in standard form. 

Proof: We first prove that Pr op er tv [T] of Definition |6 . 1 1 holds . Considering the capacity a{vn) 
of the rightmost odd edge e = {2m + 2n — 2,2m + 2n — 1} shows that F cannot contain 
both the long high profit task for Vn and the short high profit task for Vn that uses this edge. 
However, since F fully uses the capacity of e (Proposition 16.51) . exactly one of these cases is 
true. Combined with Proposition 16.61 this shows that F either contains the long high profit 
task for Vn, or all short high profit tasks. Since F is optimal, the low-profit task for Vn is 
included if and only if the long high profit task is included. This shows that Property [1] 
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of Definition 16.11 holds for Tlien tlie same argument can be applied using the odd edge 
{2m + 2ra — 4, 2m + 2n — 3} to the left of e, which has capacity a{vn) + a{vn-i), to prove this 
property for Vn-i- Continuing this way proves Property [T] of Definition 16.11 for all Vi. 

Next, we will prove that Property [2] of Definition 16.11 also holds for the optimal solution 
F. We will prove that the vertices of G for which F contains a long high-profit task form an 
independent set in G. More precisely, for any edge {vi,Vj} G E{G), F does not contain both 
the long high-profit task for Vi and the long high-profit task for vj. Let = {vi,Vj} (with 
k £ {1, . . . , m}). According to our construction this edge is associated with an even left edge 
e = {2k — 1, 2k} G E[P) of ufpp(G). The capacity of e is one unit smaller than the capacity 
of its adjacent odd edge e' = {2k — 2,2k — 1}, and yet the capacity of e' is fully used by F 
(Proposition 16. 5( ). The only tasks using e' but not e are a short task for Vi, and a short task 
for Vj. So for at least one of Vi and vj, F includes this short task. Since we already showed 
in the first part of the proof that F cannot contain both a long high-profit task and a short 
high-profit task for the same vertex Vk, it follows that for at least one of Vi and Vj, F does 
not contain the long task. Hence Property [2] of Definition 16.11 also holds for F, and thus F is 
in standard form. □ 

Combining the above lemma again with a calculation of the profit of a solution in standard 
form, this shows that any optimal solution to UFPP(G) corresponds to an independent set of 
G: 

Lemma 6.8 //UFPP(G) admits a solution of profit at least x + X]r=i ct{vi)n{m + i), then G 
has an independent set of size at least x. 

Proof: Consider an optimal solution F, so w{F) > x -\- Yll^=i Oi{vi)n{m + i). By Lemma |6.7|, 
F is in standard form. Let / be the set of vertices Vi of G for which F contains the long 
high-profit task. Since F is in standard form, / is an independent set of G. Proposition 16.31 
shows that w{F) = \I\ + '^i^iOi{vi)n{m + i). It follows that |/| > x, so / is the desired 
independent set. □ 

Our construction uses only polynomially bounded numbers (profits, demands and capac- 
ities), so even if the numbers are encoded in unary, this is a polynomial transformation. For 
this it is also essential that a 3-coloring a of G can be found in polynomial time [33j . Lem- 
mas 16.41 and 16.81 show that G has an independent set of size at least x if and only if P, F 
admits a solution of profit at least x + X^iLi 0({vi)n{m -\- i). Since the problem of finding a 
maximum independent set is NP-hard when restricted to graphs of maximum degree 3 |28j . 
this proves that UFPP is strongly NP-hard. The problem obviously lies in NP, which yields: 

Theorem 6.9 UFPP is strongly NP-complete when restricted to instances with demands in 
{1,2,3}. 

In fact, with a small addition we can show that the problem remains NP-complete when 
all edge capacities are equal. Let Um be the maximum capacity used in the instance UFPp(G) 
constructed above, and let X = n -|- a{vi)n{m + i) be an upper bound on the profit of 
any solution (Lemma 16. Sp . Note that X is bounded by a polynomial in n. For every edge 
e = {/c. A; -|- 1} of P with Ug < u^, we can increase the capacity to u^, and introduce 
dummy tasks, from k to k-\-l with demand 1 and profit X. For a polynomial transformation 
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we must ensure that the number of dummy tasks is polynomiahy bounded. Given an instance 
with demands in {1, 2, 3}, a maximum capacity of 3n suffices, so w.l.o.g. Ue < 3n holds for ah 
capacities. Hence, on each edge at most 0{n) dummy tasks are introduced. Since the number 
of edges is bounded by 0(n), at most O(n^) dummy tasks with profit X and demand 1 need 
to be added. In an optimal solution, clearly all of the dummy tasks are selected, so this yields 
an equivalent instance. Therefore the previous theorem can be strengthened to: 

Theorem 6.10 UFPP is strongly NP-complete when restricted to instances with demands in 
{1,2,3}, where all edges have equal capacities. 

7 Discussion 

In this paper, we presented the first constant factor approximation algorithm for the unsplit- 
table flow problem on a path. From a technical point of view, our key contributions are 
the framework employing the (a, /3)-approximation algorithms for the small tasks, and the 
geometrically inspired dynamic program for the large tasks. In particular, for the latter we 
establish a connection between UFPP and Maximum Weight Independent Set of Rectangles, 
two problems which at first glance might seem totally unrelated. It may be interesting to 
study whether our new techniques for finding maximum non-overlapping sets of rectangles 
can be applied to other geometric (packing) problems. 

A challenging open question is whether a PTAS exists for UFPP, or whether the problem 
is APX-hard. Note that even in the much simpler case of uniform capacities (RAP), this 
question is open; the best known approximation algorithm has an approximation ratio of 
2 + e |14j . It seems that for a PTAS one would have to use an approach that deals with all 
types of tasks (that is, small and large) in a combined manner, rather than using a best-of- 
two analysis like in our algorithm, or in other algorithms such as the (2 -|- e)-approximation 
algorithm for UFPP-NBA [20]. 

We proved strong NP-hardness for UFPP with uniform capacities (RAP), even if the 
input is restricted to demands in {1, 2, 3}. Due to this result it would be interesting to know 
whether the special case of demands in {1,2} is also NP-hard, or admits a polynomial time 
algorithm. Note that if the demands are uniform, the problem can be solved in polynomial 
time with a small extension of the algorithm by Arkin and Silverberg [4] . 

Furthermore, our research raises the question what generalizations of UFPP also admit a 
constant factor approximation. One might consider generalizations such as those that have 
been considered for RAP and UFPP-NBA, either in a scheduling context or in a network flow 
context (see Section [l.2p . In particular, an interesting open problem is whether the unsplit- 
table flow problem on trees (without the NBA) also admits a constant factor approximation 
algorithm. 
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